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ABSTRACT 

We explore whether stellar tidal streams can provide information on the secular, cosmological evo¬ 
lution of the Milky Way’s gravitational potential and on the presence of subhalos. We carry out 
long-term (At ~ thubbie) N-body simulations of disrupting satellite galaxies in a semi-analytic Galaxy 
potential where the dark matter halo and the subhalos evolve according to a ACDM cosmogony. All 
simulations are constrained to end up with the same position and velocity at present. Our simulations 
account for: (i) the secular evolution of the host halo’s mass, size and shape, (ii) the presence of 
subhalos and (iii) dynamical friction. 

We find that tidal stream particles respond adiabatically to the Galaxy growth so that, at present, the 
energy and angular momentum distribution is exclusively determined by the present Galaxy potential. 
In other words, all present-day observables can only constrain the present mass distribution of the 
Galaxy independent of its past evolution. We also show that, if the full phase-space distribution of 
a tidal stream is available, we can accurately determine (i) the present Galaxy’s shape and (ii) the 
amount of mass loss from the stream’s progenitor, even if this evolution spanned a cosmologically 
significant epoch. 

Subject headings: stellar dynamics - methods: N-body simulations- methods: semi-analytical - 
galaxies: kinematics and dynamics - galaxies: halos - galaxies: dwarfs 


1. INTRODUGTION 

In the last decade stellar streams in and around the 
Milky Way, which are possible debris from the disrup¬ 
tion of satellite galaxies during the hierarchical assembly 
of our Galaxy, have become an active topic of investiga¬ 
tion for several reasons. Firstly, large scale GGD surveys 
have provided unprecedented evidence of accretion and 
tidal disruption of dwarf galaxies around large spirals in 
the Local Group (Milky Way: see Majewski 2004; M31: 
Ibata et al. 2002 and beyond: Pohlen et al. 2004). Sec¬ 
ondly, because tidal streams provide strong constraints 
on the potential of host galaxies, it is possible to estimate 
the shape of dark matter halos on large scales in contrast 
to traditional tracers, such as HI or stellar kinematics 
(see Sackett et al. 1999 for a review), which provide esti¬ 
mates on relatively small scales. The expected shape of 
dark matter halos depends on the nature of dark matter 
particles (see, for example, Dubinsky & Galberg 1991, 
Yoshida et al. 2000 and Dave et al. 2001 for shape esti¬ 
mates for cold, self-interacting and hot dark matter mod¬ 
els, respectively), tidal streams represent a useful tool to 
discriminate between different paradigms. In addition, 
tidal stream properties also depend on the mass and in¬ 
ternal structure of its progenitor (e.g Johnston, Sackett & 
Bullock 2001 and Law, Johnston & Majewski 2004, here¬ 
inafter LJM), which ultimately results in a complemen¬ 
tary way of estimating the progenitor mass and there¬ 
fore its mass-to-light ratio. Finally, tidal streams can be 
used to determine the position of the progenitor if this 
has not previously been detected (e.g. Font et al. 2004, 
Peharrubia et al. 2005). 

The formation of tidal streams is, conceptually, a sim- 
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pie process: along its orbit, a stellar system crosses re¬ 
gions where the tidal force of its host galaxy supplies 
kinetic energy to the initially bound particles. If the 
energy gain is large enough, particles can become un¬ 
bound and escape from the host system, forming two 
sub-systems which are kinematically well differentiated: 
the leading and the trailing tails which, as their names 
indicate, precede and follow, respectively, the progenitor 
in its orbit. The orbital evolution of stripped particles is 
initially similar to that of the progenitor (e.g. Lynden- 
Bell & Lynden-Bell 1995, Johnston 1998) although, as 
they evolve in the host galaxy potential, orbits diverge 
from each other monotonically with time. However, even 
after a large number of orbital periods, tidal stream par¬ 
ticles and the progenitor system reside in well-defined 
regions of the constant of motion space (Helmi & White 
1999), therefore attesting a common origin. 

The complexity of the stream formation and evolution 
forces the use of N-body calculations in most cases. The 
existing work can be divided into a) live N-body simula¬ 
tions, where the host galaxy is formed by a given num¬ 
ber of particles initially in equilibrium and b) simulations 
where the host galaxy is represented by a non-responsive 
potential. Whereas the former takes into account the 
host galaxy’s response to the satellite, the later neglects 
this in order to save computational resources for exten¬ 
sive orbit surveys. 

However, none of the N-body simulations of tidal stream 
formation and evolution to date has accounted for the 
overall build-up of the host galaxy during the 1-10 Gyr 
that the stream formation may take. Yet, in the com¬ 
monly accepted hierarchical scenario, host galaxies ex¬ 
perience large changes in mass, size and shape during 
their history and hence a tidal stream evolving in an 
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unchanging host galaxy can only approximate recent 
epochs. The main goal of this contribution is to address 
the effect that the secular evolution of the host galaxy 
induces on the formation, evolution and interpretation 
of tidal streams. Specifically, we will explore whether 
the present-day structure of an extensive tidal stream 
can constrain the past history of the host’s gravitational 
potential. 

In addition to the secular overall mass and size growth 
of the host’s halo, we will also examine the influence of 
dark matter substructures on tidal streams (Ibata et al. 
2002 and Johnston, Spergel & Haydn 2002). CDM cos¬ 
mology predicts a large number of sub-structures in a 
galaxy-sized halo, many more than the number of ob¬ 
served dwarf galaxies (Klypin et al. 1999, Moore et al. 
1999). Recently, it has been proposed that the process of 
re-ionization in the universe would lead to a significant 
decrease in the number of “visible” sub-structures (Bul¬ 
lock, Kravtsov & Weinberg 2000, Benson et al. 2002, 
Somerville 2002, Tully et al. 2002), while keeping the 
total number of sub-structures unaltered. This process 
would establish a minimum mass (corresponding to a 
minimum circular velocity of ~ 30 km/s) above which 
gravitationally bound systems would be able to retain 
baryonic matter and, thus, to form stars. 

The properties of kinematically cold tidal stream are 
strongly sensitive to the lumpiness of the galaxy poten¬ 
tial (Ibata et al. 2002, Johnston, Spergel & Haydn 2002) 
as repeated encounters with dark matter sub-structures 
alter the energy and angular momentum distribution of 
tidal stream particles, leading to hotter, broadly dis¬ 
persed streams. Tidal streams appear to be a unique lab¬ 
oratory to determine the presence of dark matter clumps 
in galaxy halos. Firstly, because the effects of those 
clumps on tidal streams are purely gravitational (and so, 
independent of whether sub-structures enclose baryons or 
not) and, secondly, because tidal streams can be detected 
on large scales and can be as old as the host galaxy, and 
may therefore provide information on the number and 
spatial distribution of bound sub-structures at different 
epochs. Yet, the calculations mentioned above do not 
take into account the evolution of the spatial distribution 
nor the mass loss of substructures, which may weaken the 
influence of substructures on the tidal stream evolution. 

In this paper we focus on tidal streams in the Milky 
Way simply because only for our Galaxy can streams be 
resolved into stars and accurate phase-space information 
gathered. 

The remainder of this paper is arranged as follows. 
In (21 we describe our models for the Milky Way poten¬ 
tial and for the satellite galaxies that are disrupted to 
form tidal streams; details our N-body code used for 
evolving satellite orbits. SI descri bes the specific set of 
orbits we explored in this work. S31 describes how we 
address the problem of the unknown mass loss history of 
the satellite. In SI we explore the behavior of the satel¬ 
lite galaxy’s orbit while in SI we explore the properties 
of the associated tidal streams. In (Qwe explore whether 
evolution of the halo shape can influence the properties 
of tidal streams. SI examines the influence of dark mat¬ 
ter substructures of these tidal streams. Finally, in 
present our conclusions. 


2. MODELS FOR THE MILKY WAY AND THE SATELLITE 
GALAXIES 


2.1. Milky Way potential 

The host galaxy system is described by a time- 
dependent gravitational potential. Our Galaxy model 
consists of a Miyamoto-Nagai (1975) disc, a Hernquist 
(1990) bulge and a Navarro, Frenk & White (1995, 1996, 
1997) dark matter halo (hereafter, NFW halo). 

The gravitational potential of each of those compo¬ 
nents in cylindrical coordinates is: 
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where <I>d(r) and d)f,(r) are time-independent, but the 
halo potential $/t(r, t) evolves as described below. There, 
rvir,rs are the virial and scale radii, respectively, Mh = 
AF/i(rvir) and = R'^ + The potential of an 
axi-symmetric NFW halo was calculated from Chan¬ 
drasekhar (1960) using elliptic coordinates 
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where q is the axis-ratio of iso-density surfaces. In this 
contribution, we shall denote oblate and prolate halos as 
those with q < 1 and g > 1, respectively. For spherical 
halos Eq. © reduces to 
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Following Johnston et al. (1999) we fix the disc and 
bulge parameters as M^i = 1.0 x IO^^Mq, Mf, = 3.4 x 
IO^^Mq, a = 6.5kpc, b = 0.26 kpc and c = 0.7kpc. The 
Milky Way halo parameters at z = 0 were taken from 
Klypin, Zhao & Somerville (2002) being Mh = 1.0 x 
IO^^Mq, Tvir = 258kpc, Ts = 21.5kpc, which leads to a 
concentration at the present epoch of c = rvir/^s = 12. 

Additionally, we have defined a characteristic velocity, 
energy and angular momentum throughout this contri¬ 
bution, which have the values Uch = 262 km/s, Ech = 
~ 68600 (km/s)^ and Lch = 4000 kpc km/s. 

Our selection of the different Milky Way components 
leads to a Galaxy potential that is aspherical everywhere. 
Note that for the gravitational potential at the present 
epoch, the only parameter that is varied in this work is 
the halo’s density axis-ratio (g). In Fig. ^ we plot the 
flattening of the total Galaxy potential as a function of 
distance and halo axis-ratio. One can clearly distinguish 
three regions: (i) the innermost region with g$ ~ 1 dom¬ 
inated by the spherical bulge component, (ii) an inter¬ 
mediate region, 3 < r < 15 kpc, where the disk potential 
dominates, so that the Galaxy potential has an oblate 

















3 



r (kpc) 

Fig. 1.— Flattening of the total Galaxy potential as a func¬ 
tion of Galactocentric distance for different halo axis-ratios. The 
horizontal dotted line indicates t = 1. 


shape < 1) and (iii) the outer regions, r > 15 kpc, 
where the halo potential dominates. 

Some caveats must be kept in mind when using this 
simplistic Milky Way model: 

Firstly, we have not included the evolution of the bary- 
onic components. Observations at high redshift have 
shown that disc galaxies evolve in mass and size (e.g. 
Trujillo & Pohlen 2005, Barden et al. 2005), which may 
induce effects similar to those of the varying halo po¬ 
tential on tidal streams. Moreover, the disc and bulge 
formation also changes the dark matter distribution, 
mostly in the inner halo region, as has been shown by 
several authors (e.g., Dubinsky 1994, Klypin, Zhao & 
Somerville 2002). However, the self-consistent evolution 
of our Galaxy and its effects on tidal streams would re¬ 
quire the use of large N-body cosmological simulations, 
beyond the scope of this paper. Here, we shall merely 
examine tidal streams in a Milky Way-like galaxy, where 
only the dark matter component evolves. 

Secondly, this evolution of our Milky Way model ignores 
any feed-back processes that satellite galaxies may induce 
on the host system. 


2.2. Average evolution of a Milky Way-sized halo 


In cold dark matter cosmogonies, halos grow hierarchi¬ 
cally in mass and size owing to mergers of less massive 
systems. The process is stochastic in detail and is de¬ 
scribed by the merger tree of a halo. However, the av¬ 
erage properties that determine the halo density profile 
at a given epoch have been studied exhaustively, both 
analytically and numerically, showing a smooth change 
with redshift or epoch. In this framework it is possible to 
trace back the evolution of a halo from its present density 
distribution. 

For the mean mass evolution Wechsler et al. (2002) 
find 
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where a = (1 -I- z) Uc is the ’formation epoch’ of the 


halo and S = d log M/i/d log o is the log of the mass ac¬ 
cretion rate at a = Oc. S acts simply as a normalization 
constant in the fit of Mh', since its value can be arbitrar¬ 
ily chosen, we therefore fix S' = 2 as in Wechsler et al. 
( 2002 ). 

These authors also show that the average concentration 
of halos can be fit at any given time by the expression 
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For a halo concentration of c = 12 at z = 0, for the Milky 
Way this fixes the formation epoch to Oc — 0.34 (using 
S=2). 

Bullock et al. (2001) have shown that the mean evolu¬ 
tion of the virial radius can be expressed as 
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“ 1 + Z 1.0 X lOllMg/r-l lv^^oAvir(2) J ’ 

. 

where Avir is the virial over-density, which can be written 
as (Bryan & Norman 1998) 
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and n is the mass density of the Universe 
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We have assumed for this study a ACDM Universe by 
fixing Qq = 0.3, Ua = 0.7 and h = 0.7. 

Eq. 0 and 0 determine the evolution of the 
halo potential $?t(r,t) shown in Eq. 0. 


2.3. Statistical realizations of the Milky Way halo and 
the merger tree 

The stochastic growth of dark matter halos through 
mergers with smaller halos implies a tight correlation be¬ 
tween the properties of the host halos and their merger 
trees. 

In addition to the overall (smooth) mass growth out¬ 
lined in Section E21 we also analyze the effects of dark 
matter sub-structures on tidal stream evolution by in¬ 
cluding an additional force component induced by these 
systems. We generate initial conditions for our cal¬ 
culations by constructing merger trees for an ensem¬ 
ble of overall identical dark matter halos identified at 
z = 0. Specifically, we take at z = 0 dark mat¬ 
ter halos of mass 1.0 x IO^^Mq, rvir = 258 kpc and 
c = 12. Following the method of Cole et al. (2000), 
we construct a statistical realization of the merging his¬ 
tory of these halos back to high redshift, assuming a 
ACDM cosmology with Dq = 0.3, Da = 0.7, us = 0.9, 
h{= i7o/100km s“^ Mpc“^) = 0.7 and F = 0.21, and an 
inflationary initial power spectrum (i.e. P(k) oc k). At 
each point in time, we identify the most massive progen¬ 
itor of the z = 0 halo and classify this as the “host” halo 
at that redshift. We then identify halos which are about 
to merge with the host at each time and label these as 
new satellites. 

Based on such merger trees, we record the virial mass, 
virial velocity and concentration of the host halo (as¬ 
sumed to have an NFW density profile) at each epoch. 
We also record the same information for each satellite 
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Fig. 2.— Evolution of the dark matter properties for the host halo 
in our simulations. Upper-left panel: Evolution of the halo mass 
enclosed within the virial radius. Upper-right panel: Evolution of 
the halo virial radius. Bottom-right panel: Scale-radius evolution. 
Bottom-left paneLKalo circular velocity evolution (see text for the 
definition of Vc as a function of Mh-,ryj\^^rs). Thin lines show 
nine statistical realizations of the merger tree of a Milky Way-like 
galaxy. Thick lines show the mean evolution of the Galaxy halo as 
outlined in Section 


along with the epoch of each merging event. The initial 
orbit for each satellite is set by choosing an initial po¬ 
sition at random on the surface of a sphere with radius 
equal to of the host halo, while the initial velocity 

is chosen from the distribution found by Benson (2004) 
for 0 = 0 in ACDM cosmologies^. 

In this way we generate both the formation history of 
the host halo and the properties of the infalling popula¬ 
tion of satellites. These are used as inputs to our model 
for the evolution of substructures (Penarrubia & Benson 
2005). The outcome is nine different merger trees of a 
Milky Way-like galaxy for which, at each time-step, we 
have the position, velocity, mass, virial and scale radii of 
all subhalos with masses above IO^Mq. 

2.4. Satellite models 

The actual disruption of satellites is modeled by ‘live’ 
(i.e. self-consistent, self-gravitating) N-body realizations 
of satellites in an evolving, but non-responsive, halo po¬ 
tential. We also include perturbations from the ensemble 
of dark matter subhalos. 

We have realized our numerical models of satellite 
galaxies with King profiles (King 1966), which approxi¬ 
mate the light profile of the Galaxy’s low-mass satellites 
(e.g. Irwin & Hatzidimitriou 1995). The high mass-to- 
light ratios measured in dwarf galaxies (such as Ursa Mi¬ 
nor and Draco, see Mateo 1998, Kleyna et al. 2004, Lo- 
cas, Mamon & Prada 2005, Munoz et al. 2005) indicate 


that these systems are dark matter dominated , but with 
a dark matter profile that is not well determined. 

As Zhao (2004) showed, the mass and orbital evolution 
are fairly sensitive to the specific mass profile we assume 
for our satellites, an aspect that must be kept in mind, 
but it is not pursued further here. 

For the initial satellite structure we assume a King 
profile with a dimensionless central potential Wq = 4, or 
concentration parameter c = logiQ^rt/rk) — 0.84, where 
Tk and rt are the King and tidal radii, respectively. Our 
satellite models have 10® particles. The initial mass and 
tidal radius are chosen such that the mass at the present 
epoch is Ms(t = tj) = 5 x 1O®M0 in all of our calcu¬ 
lations. Details of how to make this choice are given in 
Section^ 


3. MODELING THE DISRUPTION OF N-BODY SATELLITE 
GALAXIES 

The force acting on each satellite particle in our simu¬ 
lations has four components: 1) the gravity of all other 
satellite particles; 2) the force from the smooth, non- 
responsive Milky Way potential; 3) dynamical friction 
on the still bound satellite portion and 4) the time- 
dependent forces from subhalos. 

To calculate the force from self-gravitating satellite 
particles we use Superbox, which is a highly efficient 
particle mesh-algorithm based on a leap-frog scheme (for 
details see Fellhauer et al. 2000). 

The acceleration from the Milky Way gravitational po¬ 
tential is calculated from Eq. o, m and © and in¬ 
cluded in the N-body code as an external force. 

The code explicitly includes the dynamical friction 
force term resulting from the interaction of the satel¬ 
lite galaxy with the smooth dark matter component. In 
axi-symmetric self-gravitating halos, the presence of a 
satellite galaxy traveling through the system induces a 
density wake behind the satellite, which exerts a drag 
force on the satellite particles. Following Binney (1977) 
this can be approximated by 
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where i = x,y and {aa^az) is the velocity dispersion el¬ 
lipsoid in cylindrical coordinates with constant ellipticity 
= 1 — (cTz/tTfi)^. Here Br and Bz are given by 
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with (vRjVz) is the satellite velocity in this frame and 
In A the Coulomb logarithm. 

The NFW density profile in elliptical coordinates can 
be written as 


^ Benson (2004) finds some evidence for evolution of the velocity 
distribution to higher redshifts, and for dependence on the masses 
of the merging halos. These dependencies are poorly quantified, 
and so we ignore them here. 
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where we have that m^(0) = {R^ + z^lq^)lrl from 

Eq. 0). 

At each point in time, the satellite mass Mg is defined 
as the sum of still-bound particles (see Section Oil . 

Peharrubia, Just & Kroupa (2004) have checked these 
equations against several satellite orbits in self-consistent 
axi-symmetric systems by keeping the Coulomb loga¬ 
rithm as a free parameter. Best fits were found for 
InA = 2.1, almost independent on the satellite’s mass, 
orbital eccentricity and halo axis-ratio. 

In order to implement Eq. CD, we have assumed that 
only bound particles feel dynamical friction, whereas un¬ 
bound (stripped) particles do not. In this way we are 
neglecting the effects that the self-gravity of the den¬ 
sity wake may induce on the orbit of the escaping par¬ 
ticles. This appears a sensible approximation as most 
of the unbound particles remain close to the main sys¬ 
tem (and hence the density wake) for only a short time 
and, furthermore, that /df ~ [10“ , 10“^]d$r/(ir (Just & 
Peharrubia 2005). 

The last force component results from the presence of 
dark matter clumps. Using the semi-analytic algorithm 
presented in Peharrubia & Benson (2005), we calculate 
and store the position, velocity and internal properties 
of dark matter sub-structures at each time-step. This 
semi-analytic code has been proved to reproduce fairly 
well the mass and spatial distributions of dark matter 
sub-structures obtained from large cosmological N-body 
simulations. In each step of the satellite simulation, the 
force on the from all subhalos is 


f. , _ __ 

■*■2, sub — — \ / 9 , I l9\‘:i/9 

^ (e^ -h |rj - r,|2)3/2 


(13) 


where Cj is the softening parameter of the substruc¬ 
ture which accounts for its finite extent. We assume that 
Cj = r/jj, where rhj is the half-mass radius. 

We must remark that, in contrast to previous studies of 
the effects of halo lumpiness on tidal streams (Ibata et al. 
2002, Johnston, Spergel & Haydn 2002), our approach 
accounts for the evolution of mass, spatial distribution 
and size of dark matter sub-halos. 

After calculating all the force terms, Superbox solves 
the equation of motion of each satellite particle 


dt"^ 


— V ($s -|- <I>d -|- d>i, -f ^h)i + + fii,sub (14) 


where is the self-gravitational potential of the satel¬ 
lite galaxy calculated by superbox, ^ = 0 for unbound 
particles and C = 1 for bound ones. Superbox uses a 
leap-frog scheme with a constant time-step of At = 0.65 
Myr, which is about 1/lOOth the dynamical time of our 
most massive satellite model. We have three resolution 
zones, each with 64^ grid-cells: (i) The inner grid cov¬ 
ers out to 10 King radii, providing a resolution of r^/G. 
(ii) The middle grid covers 20 King radii r^/S. (iii) The 
outermost grid extends to 348 kpc and contains the local 
universe, at a resolution of 11.6 kpc. At any time-step, 
all grids are centered at the density maximum of the 
satellite model. 

This choice of grids provides maximal resolution within 
the satellite itself, where the self-gravity of the satel¬ 
lite dominates. Most unbound particles eventually reach 
r > 20rfe from the satellite center, i.e, within the lowest- 
resolution zone. This, however, induces negligible effects. 


for the orbit of stripped particles is mainly determined 
by the semi-analytical Galaxy potential (and so is insen¬ 
sitive to spatial-resolution problems). 

4. CALCULATIONS 

With the modeling tools at hand we now explore what 
can be learnt about the dynamical history of a disrupt¬ 
ing satellite and its tidal tail on the basis of present-day 
observables. One additional open point addressed in this 
Section is how to reconstruct the mass-loss history and 
the orbit of a satellite galaxy using an iterative method. 

4.1. Present-day parameters 

We have focused our study on a Sagittarius-like galaxy 
because, at the present day, it is the system in the 
Milky Way for which the most comprehensive observa¬ 
tional constraints exist. Eollowing LJM (and references 
therein) we adopt the following remnant properties at 
t = tf = 14Gyr: 

• Mg{tf) = 5 X 10®M©; 

• r{tf) = Tperi = 15kpc; 

• Vte,n{tf) = 320 km/s; 

• Orbital inclination i = 45°. 

With the present mass, position and velocity of a satel¬ 
lite galaxy fixed, we explore the physical processes that 
may have influenced the orbit of the satellite galaxy in 
the past: 

1) Axis-ratio of the host halo. 

2) Dynamical friction. 

3) Smooth evolution of the Galaxy’s potential. 

4) The presence of a large number of dark matter sub¬ 
structures. 

To analyze the effects of various factors on the satellite 
orbit and stream evolution we have conducted several 
N-body simulations. Tabled lists the properties of the 
satellite and host halo used in each simulation. 

4.2. Iterative reconstruction of the satellite orbit 

Each model has been evolved in the following way: 

(i) The present time is denoted as tf = 14 Gyr. Start¬ 
ing from the present, we integrate models backwards 
in time to to = 4 Gyr ^ solving Eq. El with 
^ = — 1 and fsub = 0 for a single particle. The goal 
is to determine r, v at t = to for a given mass loss history. 

(ii) Using the results of (i), we place each N-body 
satellite realization at r(to) with velocity —v(to). Subse¬ 
quently, we evolve the system forwards in time using the 
N-body algorithm selecting those realizations which 
produce N-body satellites that match our constraints at 
t = tf . 

This method assures that all models match the present 
constraints independently of their past evolution. How¬ 
ever, point (i) cannot be trivially solved since the solution 
of the equation of motion (Eq. [T^ is coupled to that of 

^ In order to fix to we have applied the results of Zentner & 
Bullock (2003) and Peharrubia & Benson (2005). These authors 
show that, as a result of tidal disruption, the maximum accretion 
time of sub-structures with M{tf) > 10^ Mq is — to — 10 Gyr. 
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TABLE 1 
Models 


Name 

MIHIQI 

M1H1Q2 

M1H1Q3 

M1H2Q1 

M1H2Q2 

M1H2Q3 

M2H1Q1 

M2H1Q2 

M2H1Q3 

M2H2Q1 

M2H2Q2 

M2H2Q3 

Ms{to)/Ms(tf) (1) 

2 

2 

2 

2 

2 

2 

10 

10 

10 

10 

10 

10 

Tk (kpc) (2) 

0.55 

0.55 

0.55 

0.58 

0.58 

0.58 

1.20 

1.20 

1.20 

1.30 

1.30 

1.30 

<? 

0.8 

1.0 

1.2 

0.8 

1.0 

1.2 

0.8 

1.0 

1.2 

0.8 

1.0 

1.2 

Halo (3) 

static 

static 

static 

evolving 

evolving 

evolving 

static 

static 

static 

evolving 

evolving 

evolving 


(1) We fix io = 4 Gyr and tf = tnubble = 14 Gyr. All satellite models have a final time mass of M(ty) = 5 X 10®Mq. 

(2) King radius at f = to- Fol Wo/cr'^ = 4 the tidal radius is rt = 6.9rfc. 

(3) Static halos have M(,(t) = M^(tj-), r’s(t) = Tsitf) and 


the mass evolution through the dynamical friction term 
(which scales as pg{r)Ms{t), where pg is the Galaxy den¬ 
sity profile). In practice, that means that the satellite’s 
orbit and mass must be simultaneously calculated hack- 
wards in time. 


4.3. Modeling satellite mass loss 

In order to disentangle the past mass and orbit evolu¬ 
tion we have applied the following scheme: 

The fist step is to fix the total amount of mass loss. Ob¬ 
servational data indicate that satellite galaxies appear 
dark matter dominated (e.g. Mateo 1998), which pre¬ 
cludes a reconstruction of the total mass loss history 
from the present-day distribution of stars within a satel¬ 
lite galaxy. Therefore, we treat the total amount of mass 
loss as a free parameter. According to Peharrubia & 
Benson (2005): (i) The present-day population of sub¬ 
halos have lost on average 50% of their mass since they 
were accreted and (ii) mass loss events where subhalos 
lose more than 90% of their accretion mass are extremely 
rare. To analyze typical and extreme cases we carry out 
simulations where Ms{to) = 2Ms{tf) = 10 ^Mq (models 
Ml) and Ms{to) = 10Ms(tf) = 5 x IO^Mq (models M2, 
see Table ^|. 

The second step is to determine the mass loss history, 
i.e, Ms{t). Here we use the results of Zhao (2004), who 
showed that, for a given satellite density profile, the mass 
evolution of satellite galaxies can be described by a set 
of empirical functions. In particular we use 


Ms{t) 

Msitf) 


Ms{to) 

Msitf) 



-I- < exp 




tf ) 


(15) 


where the power p depends on the mass loss fraction 
Msitf)/Msito). The time interval is t G [to,tf] with 

tf<tQ. 

Lastly, Eq. m still requires to fix the value p for each 
simulation. Since orbit and mass are uniquely related 
once Msito)/Msitf) has been chosen, we use the follow¬ 
ing iterative scheme to determine p and solve Eq. m 
for each satellite model: 

0) Selecting the starting p-value: Assuming Ms(t) = 
const., we integrate the orbit backwards in time. Sub¬ 
sequently, we run the N-body orbit forwards in time in 
order to calculate the number of bound particles as a 
function of time and to fit an initial, tentative value of p. 



4 6 a 10 12 14 


1 (Gyr) 

Fig. 3. — Mass evolution of our satellite models. Upper and 
lower panels show the mass loss of satellites with initial masses 
Msito) = 5 X IO^Mq (family model M2) and Msito) = IO^Mq 
( family model Ml), respectively. All models have the same mass 
at present (Ms[t/] = 5 x 10 ®Mq). Thick lines show the curves 
obtained from Eq. 113. Dotted lines show the fraction of bound 
particles obtained from the self-consistent satellite potential cal¬ 
culated by SUPERBOX for each satellite of the same family model. 
Note that the mass loss curves only depend on Msito) / Msitf) and 
not on the exact Galaxy potential. 


\)Orhit calculation: With the tentative 75 -value we inte¬ 
grate Eq. d backwards in time with Ms(t) given by 

Eq. (IT^ . 

2) Fixing the initial satellite size and the total mass loss: 
We have used the fact that the initial satellite’s tidal ra¬ 
dius and its mass loss rate are correlated. Using the mass 
evolution and the orbit obtained in (i), we run several N- 
body simulations varying the initial rt (see with a 
fixed Mg (to) and select those models that end up with 
Msitf) ~ 5 X I0®Mq. 

3) Fine-tuning: Finally, we fine-tune the value of p to 
obtain a better match the the mass-loss histories derived 
from the N-body simulations. 

This process repeats until the N-body and the empiri¬ 
cal Ms(t) match^, typically 2-3 iterations. 

® Note that Ms(tf) does not enter into Eq. and so the 

present mass depends only on the value of p. Since p is found 
from the best fit to the mass curve within [to,ty], Eq. 1151 may, 


























7 


In Tabled we have listed the properties of our satel¬ 
lite models and in Fig. 13 we show the mass evolution of 
satellites with Ms{to)/Ms{tf) = 10 (upper panel, family 
models M2) and Ms(to)/Ms{tf) = 2 (lower panel, family 
models Ml). Thick lines show the empirical evolution ob¬ 
tained from Eq. 113 with p = 1 (lower panel) and p = 4 
(upper panel). Dotted lines show the amount of bound 
mass calculated from the N-body models as a function of 
time: The number of bound particles at any time-step is 
determined from the self-consistent potential calculated 
by SUPERBOX. We label a particle as unbound whenever 
E = 5 (v — Vcm)^ + ‘I’s >0, where Vcm is the center- 
of-mass velocity of the bound particles, v is the particle 
velocity in the Galaxy frame and is the potential in¬ 
duced by the satellite’s self-gravity. 

Interestingly, we find that the shape of the mass 
loss function only depends on the total amount of 
mass loss and not on the dark matter halo’s parame¬ 
ters, such as mass, size and density axis-ratio, so that 
once Ms{to), Ms{tf), (rs,Vs)(t/) are fixed, there exists a 
unique value rt(to) that reproduces the function Ms(t) 
shown in Eq. 11511 . 

As commented above, the exact value of rt(fo) and 
Ms{tQ)/Ms{tf) are directly correlated. As we show in Ta- 
blem satellite models in evolving halos (models H2) have 
initial tidal radii a factor 1.06 larger than those models in 
static halos (models HI), which indicates that mass loss 
is slightly reduced if halo evolution is implemented. The 
relation between the initial tidal radius and the mass loss 
fraction is more evident if comparing models Ml against 
models M2. There we see that strong mass loss events 
imply large initial tidal radii. In particular we find that 

rt,M2/n,Mi - 2-2 for = 5 at f = to- 

Although this technique is computationally expensive 
(for every simulation one must run several N-body sim¬ 
ulations in order to fix rt(to) and p), it avoids the use of 
complex semi-analytic algorithms to estimate mass evo¬ 
lution (e.g. Taylor & Babul 2001, Zentner & Bullock 
2003, Penarrubia & Benson 2005) which are considerably 
less accurate. 

5. SATELLITE ORBIT 

Here we describe the main factors that determine 
the long-term orbit evolution of the satellite up to the 
present. 

5.1. Dependence on the Halo Axis-ratio 

As already mentioned in the Introduction, the Galaxy’s 
halo axis-ratio has not yet been fully constrained. In 
practice that means that, whereas the present angular 
momentum can be directly measured from observations, 
the orbital energy has a range of possible values depend¬ 
ing on the halo’s axis-ratio (g). 

In Fig. 0] we show the value of the energy at t = tf 
as a function of q (right panel) and the resulting last 
apo-center distance {t ~ tf, left panel). As we can see, 
the apo-center is fairly sensitive to the halo axis-ratio, 
specially for g < 1 , which allows this last quantity to be 
constrained from the spatial distribution of the “young” 

for large values of Ms (to), give a slightly different value for Ms(ty) 
than desired. We have checked that this small mismatch induces 
negligible effects on the orbit integration. 


.4 0,6 0,8 1 1,2 1,4 1, 

q 




Fig. 4.— Left panel: First apo-center distance as a function of 
the halo axis-ratio (q). Right panel: Initial energy as a function of 
q. In all calculations we fix the present satellite position at r = 15 
kpc (peri-center) with a tangential velocity of titan = 320 km/s. 


tidal stream^. In Section El we shall discuss other con¬ 
straints on g from long-term orbital properties (i.e, from 
older stream pieces). 

5.2. Satellite orbit and the Mass-Loss History 

The orbit and the mass of satellite galaxies are tightly 
correlated through dynamical friction. 

In Fig. El we plot the galactocentric distance as a func¬ 
tion of time in the absence of dynamical friction (upper 
panel) and for two different mass loss histories: Models 
Ml (middle panel) and M2 (bottom panel). For clarity, 
we only show orbits in a spherical halo. Solid and dotted 
lines denote, respectively, orbits in static and evolving 
halos. 

Looking at the orbits in static halos, we find that 
dynamical friction decreases both the peri- and apo- 
center distances, which for the model Ml is approxi¬ 
mately rapo(to)/?'apo(t/) — 1-33 and for the model M2 
rapo(to)/T'apo(t/) — 2 (similar values are found for rperi)- 
Note that the orbit of the model M2 decays only 1.5 
times more than Ml, even though its mass loss is five¬ 
fold. The fact that the orbital decay is not simply 
proportional to the satellite mass can be readily shown 
by a dimensional calculation: taking into account that 
fdf oc Msp ^ Ms/r^, we have that /d/,M2//d/,Mi ^ 
Ms,M 2 /Ms,Mir^po,Mi/’’apo,M 2 - 1-48 at t = to (note that 
this value is fdf,M 2 /fdf, mi = 1 at t = tf) which indicates 
that dynamical friction rapidly decreases as the averaged 
galactocentric distance increases. 

5.3. Satellite orbit and the Evolution of the Milky Way 

Halo 

In Fig. El we also plot the orbit evolution of the models 
Ml and M2 in an evolving halo (dotted lines). Because 
the halo was smaller and less massive in the past, the 
apo and peri-center distances were larger at earlier times 
than at t csi tf, even in the absence of dynamical friction 
(upper panel). Gomparing orbits in static and evolving 

^ As an example, LJM have used the radial extension of the 
Sgr stream to perform a similar study. Unfortunately, no accurate 
proper motions of the Sgr main system are available, so that utan 
was a free parameter in their models. 
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Fig. 5.— Galactocentric distance as a function of time. In the 
upper panel we show the satellite’s orbit in absence of dynamical 
friction. In the middle and bottom panels we show, respectively, 
the orbits of the models Ml and M2. Solid and dotted lines denote 
orbits in static and evolving halos, respectively. The axis-ratio is 
fixed to q = 1. Note that all models end up at present (t = 14 
Gyr) at the same position and with the same velocity. 


halos, we find that the maximum increase in the apoc- 
entric distance (i.e. at t = Iq) is approximately 10%, 
independent of the satellite mass. 

Although, at first glance, halo evolution and dynami¬ 
cal friction appear to induce similar effects, the orbital 
properties evolve in quite different ways. In Fig. IHlwe 
plot the evolution of the energy and angular momentum 
in a static and an evolving Galaxy potential for the mod¬ 
els shown in Tabled 

The upper row shows orbits in the absence of dynamical 
friction. Without dynamical friction in a static halo, the 
energy and the perpendicular component of the angular 
momentum (L^) are constants of motion, so that a given 
orbit is represented by a point in the E — plane®. The 
different energy values for q = 0.8,1.0,1.2 correspond to 
those plotted in Fig.0 

In contrast, if the halo potential evolves, is an adi¬ 
abatic invariant (see §3.6 of Binney & Tremaine 1987) 
while the orbital energy increases (note that A < 0) as 
we integrate the orbit backwards in time, resulting in an 
increase of the peri and apo-centric distances shown in 
Fig. El (upper panel). 

If dynamical friction is switched on (middle and bottom 
rows) both, E and L^, had higher values in the past. 
Comparing curves with different halo axis-ratios, we find 
that, independent of the satellite mass, the secular shift 
of energy and angular momentum is larger in prolate ha¬ 
los (<7 > 1) than in oblate (g < 1) ones. 

Therefore, we conclude that, even without dynamical 
friction, the halo growth does not alter the momentum 
evolution, as this last quantity is an adiabatic invariant, 
although the energy variation increases by a factor ~ 1.7 
for all satellite models, independently of satellite mass 

® In practice, we see a short, horizontal line in the figure as 
energy is not perfectly conserved in our N-body scheme. 


Fig. 6.— Evolution of the satellite’s orbit in the constant-of- 
motion plane. Upper, middle and bottom raws show, respectively, 
orbits in the absence of dynamical friction and the orbits of models 
Ml and M2. Dotted, full and dashed lines show orbits in halos 
with axis-ratios of q = 0.8,1.0,1.2, respectively. The arrow in the 
bottom-right panel indicates the direction of the time evolution. By 
construction, all models have the same mass, position and velocity 
at t = t y. 


and halo axis-ratio. 

6. TIDAL STREAM PROPERTIES 

In Section El we have analyzed the orbits of satellite 
galaxies under different mass loss histories and halo po¬ 
tentials, with the only constraint that the present mass, 
velocity and position is the same in all models. In this 
Section we use N-body simulations to explore the prop¬ 
erties of the tidal streams associated with these satellite 
galaxies. The goal is to determine how much information 
on the halo potential and on the satellite mass history 
can be obtained from tidal streams. 

6 .1. The halo flattening and the orbital plane evolution 

At any instant, the angular momentum L of a particle 
orbiting in a gravitational potential corresponds to the 
normal vector of the orbital plane. Therefore, in spher¬ 
ical coordinates, the orbital plane can be described by 
two angles, the azimuthal angle (j) = ta,n~^ {Ly/L^) and 
the orbital inclination i = cos~^{Lz/\L\). If the Galactic 
potential is not spherical and i ^ 0° and i ^ 90°, we 
expect the orbital plane to nutate {di/dt ^ 0) and to 
precess {d(j)/dt ^ 0). 

As Johnston, Law & Majewski (2005) have shown, de¬ 
termining the precession rate of tidal streams provides a 
powerful technique to measure the Galactic shape since 
both effects, precession and nutation, are induced by the 
quadrupole of the total Galaxy potential. Moreover, this 
method does not require any kinematical information, 
since the normal vector can be obtained by fitting the 
3D distribution of debris onto a plane. 
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Age (Gyr) Age (Gyr) 

Fig. 7.— Precession angle {(f)) between the tidal debris and 
the parent satellite at t = tf as a function of the stream age, 
for different halo axis-ratios and initial {t = to) satellite masses. 
Filled and and open symbols denote particles in the leading and 
trailing tails, respectively, (pp is the present azimuthal angle of the 
stream’s progenitor. Solid lines show the progenitor’s precession. 
These models assume a static halo (models HI). 


In Fig.[7|we show the evolution of the azimuthal angle 
as a function of the stream age® for the models HI (static 
halo, see Tabled). We have divided the stream particles 
of the present epoch into leading and trailing tail par¬ 
ticles (filled and open symbols, respectively) for differ¬ 
ent age sub-samples. Subsequently, we calculate {(p) and 
the mean variance (this last represented by error bars). 
Additionally, we have plotted the evolution of the pro¬ 
genitor’s azimuthal angle (full lines). This Figure shows 
that: 

(i) Tidal streams do not present a single orbital plane 
but two, corresponding to the trailing and leading tails. 

(ii) As Johnston et al. (2005) finds, the precession rate 
dcj)/dt has a different sign for leading and trailing tail 
planes. In particular, leading tail planes precess as one 
would expect in an oblate galactic potential, whereas 
trailing tail planes precess in the opposite sense. 

(iii) The precession rate is sensitive to the total amount of 
mass loss. Satellite models that have lost a large amount 
of the initial mass show a smaller precession rate in the 
oldest parts of the stream in comparison with satellite 
galaxies that were originally less massive. That reflects 
the fact that, as a result of a larger dynamical friction, 
models M2 moved at larger distances than models Ml 
(see Fig. 0 ), in regions where the Galaxy potential is 
more spherical (see Fig.0|. 

(iv) Leading tails show a larger dispersion around the 

mean precession angle than trailing tails. We can also 
see that the variance (shown with error bars and defined 
as {{(j) — of the leading tails increases with the 

initial satellite mass and with the stream age. 

(v) The evolution of the satellite’s azimuthal angle is 
fairly similar to that of the leading tail if the Galaxy 

® We define “age” of a stream particle as tf —tu, where tu is the 
time when the particle became unbound and tf = lA Gyr (present). 



Fig. 8.— Width of the tidal tails (left columns) and their velocity 
dispersion (right column) a,t t = tf, measured perpendicular to the 
orbital plane as a function of age. Solid and dotted lines denote 
particles in the leading and trailing tails, respectively. Models Ml 
and M2 are shown with thin and thick lines, respectively. These 
models assume a static halo. 


halo is oblate. If the halo is prolate, the precession rate 
has opposite sign, as expected, and the trailing tail fol¬ 
lows the precession of the satellite’s orbit. 

(vi) As Johnston et al. (2005) have shown, the preces¬ 
sion of the orbital plane is highly sensitive to the halo’s 
shape. Defining A(() = {(j))i — {(j))t, where {(p) i , {p)t are the 
azimuthal angles of the leading and trailing tails, we find 
that dAp/dt ~ 40°, 25°, 12° per Gyr for q = 0.8,1.0 and 
1.2, respectively (in contrast, the progenitor’s precession 
rate was 30°, 14°, -2° per Gyr). 

6.2. The satellite mass and the stream’s width and 
velocity dispersion 

As shown by Johnston (1998), the width and the ve¬ 
locity dispersion perpendicular to the orbital plane of the 
youngest stream parts provide strong constraints on the 
mass of the progenitor. As the age increases, streams 
become wider and “hotter” (e.g. Helmi & White 1999), 
with properties reflecting the shape of the galaxy poten¬ 
tial. 

In Fig. 0 we plot the width (az) and the velocity 
dispersion (aw) in the direction perpendicular to the 
stream’s orbital plane (denoted here as the Z-axis) as 
a function of age. These quantities are defined simply as 
a^ = {(x—{x))^)^^^. We divide the sample of stream par¬ 
ticles at present into leading and trailing tails (solid and 
dotted lines, respectively) and into sub-samples of dif¬ 
ferent ages. For each sub-sample, we calculate the mean 
orbital plane in order to determine its perpendicular vec¬ 
tor. We use thin and thick lines to represent values from 
the models Ml and M2. This Figure shows that: 

(i) Young tidal streams (age < 1 Gyr) present similar 
width and velocity dispersion values, az — 2-3 kpc and 
aw — 10-15 km/s, barely dependent on the progenitor’s 
mass evolution and on the halo’s axis-ratio, in agreement 
with Johnston (1998). 

(ii) Older stream parts show a clear increase of width and 
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velocity dispersion as a function of age. That increase 
is larger in streams originating from massive progenitor 
systems and from progenitors orbiting in an oblate halo. 
Distinguishing between leading and trailing tails we ob¬ 
tain that the former are wider than the later for age > 2 
Gyr. Curiously, all models show that for 1 < age < 2 
Gyr the trailing arm is slightly wider than the leading 
one. 

(hi) The stream’s velocity dispersion has a complex 
evolution. Whereas the trailing tail shows dcrw/df ~ 
const > 0 for the whole age range, the leading tail has a 
steep rise of aw followed by a plateau at age~ 5-6 Gyr. 
The value of aw at the plateau is comparable to the 
velocity dispersion of halo stars (see e.g. Binney & Mer- 
rifield 1998), which indicates that the kinematics of stars 
that belong to the oldest stream pieces are dominated 
by the Galaxy potential, regardless of their “external” 
origin. 

As Ibata et al. (2001) showed, the thickness and veloc¬ 
ity dispersion are directly correlated with the precession 
rate of the tidal stream. By comparing Fig. [7| and El we 
can readily see that larger precession rates always induce 
tidal streams to become thicker and “hotter”. 

In the absence of proper motions, tidal streams are usu¬ 
ally identified from the background stellar population as 
spatial over-densities with a “cold” (i.e small aw) ve¬ 
locity distribution. Our results appear to indicate that, 
whereas the trailing tails can be easily detected, the old 
parts of the leading tails might be extremely hard to dis¬ 
tinguish from halo stars if the only constraints available 
are the two mentioned above. 

6.3. Streams in the E — plane. Is secular halo 
evolution detectable? 

Full observational coverage of a tidal stream is a dif¬ 
ficult task, firstly because tidal streams are spread out 
over a large range of distances and angles (therefore re¬ 
quiring large observational surveys, such as SDSS, RAVE 
or GAIA) and secondly because it is difficult to distin¬ 
guish between tidal stream and Galactic stars, especially 
for old parts of the stream. As Helmi & White (1999) 
showed, tidal streams populate well-defined regions of 
the constant-of-motion space, independent of their spa¬ 
tial distribution and age, which can enormously facilitate 
their identification. 

In Fig. El we show the present location in the 
energy-angular momentum plane of the N-body models 
M1H1Q2 (upper panel) and M2H1Q2 (bottom panel) at 
t = tf. Both models were evolved in a spherical, static 
halo. Grey dots show particles that remain bound at the 
end of the simulation, whereas black dots denote tidal 
stream particles. Tidal streams show a bimodal E and 
Lz distribution, corresponding to the leading and trail¬ 
ing tails. In comparison with the remnant system, lead¬ 
ing tails have lower E, mean values and trailing tails 
higher values. 

Satellite galaxies that were initially more massive have 
suffered a larger orbital decay through dynamical friction 
and show, therefore, a larger spread in the constant-of- 
motion plane (compare model M2 against model Ml). 
Remarkably, the area in the E — Lz plane occupied by 
bound particles only depends on the present mass and 
not on the mass loss history. That fact demonstrates 
that the mass evolution of a bound stellar system can 



Fig. 9.— E — Lz distribution of the particles in the satellite 
and the two tidal tails a.t t = tf, for model M1H1Q2 (upper panel) 
and M2H1Q2 (bottom panel). Solid and dotted lines show the 
distribution of tidal stream and bound particles, respectively. All 
distributions have been normalized to the number of satellite par¬ 
ticles (N = 10®). 


be hardly estimated exclusively from the present E-Lz 
distribution of bound particles. 

We have separated the present distribution of tidal 
stream particles into bins of different ages and measured 
the mean energy and angular momentum of the leading 
and trailing tails. In Fig. 1101 andllllwe plot those quanti¬ 
ties for all our satellite models. Thick lines also show the 
look-back evolution of the progenitor’s energy and angu¬ 
lar momentum (solid and dotted lines denote models in 
static and evolving halos, respectively). In this Figure 
we can see that: 

(i) The evolution of the progenitor’s energy and angular 
momentum induced by dynamical friction can be eas¬ 
ily traced from tidal streams if we know the time when 
stream particles were stripped out. This ultimately al¬ 
lows us to determine accurately the mass history of satel¬ 
lite galaxies with associated tidal streams. 

(ii) The present energy and angular momentum of stream 
particles are not sensitive to the adiabatic growth of the 
dark matter halo, independent of the time when stream 
particles were stripped from the progenitor. 

As shown in Fig. iini the angular momentum of the main 
system is independent of whether the halo evolves or not, 
since Lz is an adiabatic invariant, which is also applica¬ 
ble to stream (unbound) particles. In contrast. Fig. ^2 
shows that the energy of the main system followed a fairly 
different evolution in evolving and static halos (compare 
thick dotted and solid lines). However, that is not re¬ 
flected in the present energy distribution of the stream, 
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Fig. 10.— Averaged angular momentum at t as a func¬ 

tion of stream age. Error bars represent the variance around the 
mean value. Open and filled symbols denote stream particles that 
have orbited in static and evolving halos, respectively. Squares 
(circles) represent trailing (leading) particles. Solid (dotted) thick 
lines show the angular momentum of the progenitor’s orbit as a 
function of look-back time for simulations in static (evolving) ha¬ 
los. Note that particles in evolving and static halos show the same 
angular momentum distribution because Lz is an adiabatic invari¬ 
ant. 
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Fig. 11. — Average energy at t = tf as a function of stream age. 
We use the notation of Fig. ITol Note that particles in evolving 
(filled symbols) and static (open symbols) halos show at present 
the same energy distribution, even though the satellite orbital en¬ 
ergy was different in the past (compare dotted against solid lines, 
respectively). 


PARTICLE ENERGY ANGULAR MOMENTUM 



Fig. 12.— Energy (left) and angular momentum (right) distri¬ 
bution of the simulations M2H1Q2 (solid lines; static halo) and 
M2H2Q2 (dotted lines; evolving halo) at different snapshots. All 
curves are normalized to the number of stream particles at the 
given time. 


independent of its age (compare filled and open symbols). 

The reason is that the energy of the stream particles 
re-adjust at any time-step to the evolution of the halo’s 
potential. This can be clearly seen in Fig. ^1 where 
we plot the E, distributions at < = 8,10,12 and 14 
Gyr. Solid and dotted lines show the distributions from 
the models M2H1Q2 (static halo) and M2H2Q2 (evolv¬ 
ing halo). As we commented above, is insensitive to 
the time-dependence of the Galaxy potential. In con¬ 
trast, the energy distribution of stream particles varies 
as a result of the halo’s growth. 

This figure shows that, at any given time, the mean value 
of the stream energy is approximately that of the pro¬ 
genitor system (see thick lines in Fig. Illll . As a result, 
the present energy of stream particles is determined ex¬ 
clusively by the present Galaxy potential and do not de¬ 
pend on the way the Galaxy has evolved. In other words, 
tracking back the halo evolution using tidal streams is a 
degenerate problem, as we have only information at the 
present day. 

6.4. Spatial distribution of debris 

Two tidal streams with the same E-L^ distribution 
and with their progenitor systems at the same position 
will present the same spatial and kinematical properties. 
In order to show that, we plot in Fig. the spatial 
projection into the orbital plane of the present distri¬ 
bution of debris for the models M2H1Q2 (static halo, 
q = 1, left panels) and the model M2H2Q2 (evolving 
halo, q = I, right panels). We have divided each model 
into leading (grey dots) and trailing (black dots) particle 
samples. Additionally, each row shows the particle dis¬ 
tribution for different stream age intervals (given in the 
upper-right corner in gigayears). 

This figure shows interesting features: 

(i) Leading tail particles move, on average, at closer 
Galactocentric distances than trailing tail particles, as 
one would expect from their lower energy and angular 



























12 


-' ' ' ' 1 ' ' 1 1 1 ' ' ' ' 1 1 1 1 1J 
STATIC [0 , 2.S]: 

r HALO 

• ' ' 1 1 1 ' ' 1 1 i ' ' 1 1 1 ' ' ' ' - 
EVOLVING [0 , 2.5]: 

r HALO 7 

[3.5 . 5]: 

: • ; 

^ , 1 , . , , 1 , , , , 1 , , . ,' 

[2.5 , 5]: 

: # ! 

[5 , 7.5], 

L . 7.5li 

[7.5 .lOJ 

r ^ - 

[7,5 .lOj 

^ i 


-TOO D 100 

(kpc) 


-100 0 1 00 
(kpol 



Fig. 13.— Projection into the present orbital plane of the present 
distribution of tidal stream particles for the models M2H1Q2 
(static halo with q = 1, left column) and M2H2Q2 (evolving halo 
with q = 1, right column). Grey and black dots represent trailing 
and leading tail particles, respectively. We have divided the parti¬ 
cle sample into different age intervals (see upper-right corners, the 
age is given in gigayears). 


momentum. 

(ii) Trailing tails keep their coherent structure even if 
they were stripped 10 Gyr ago, whereas leading tails dis¬ 
perse within 2.5 Gyr, making their detection as spatial 
over-densities more difficult. If we compare these results 
with those shown in Fig |7| we come to the conclusion 
that the “broadening” and “heating” of tidal streams is 
more prominent in systems that suffer a large precession, 
(hi) As one expects from the results of >1(1.hi the present 
spatial distributions of debris from progenitors that have 
orbited in static and evolving halos are essentially indis¬ 
tinguishable if we force the progenitor to end up with the 
same mass, position and velocity for all our models. 

7. EFFECTS OF THE HALO SHAPE EVOLUTION 

In this contribution we have analyzed the effects that 
the growth in mass and size of the Galaxy induce on 
the stream’s properties. However, we have assumed that 
the shape of the halo remains constant in time. Accord¬ 
ing to recent cosmological simulations, that assumption 
appears to be false. Dissipationless AGDM simulations 
show that, as a result of accretion, dark matter halos be¬ 
come more spherical in time (see e.g. Allgood et al. 2005 
and references therein). For a Milky Way-like galaxy, 
these authors show that {q){z = 3) ~ 0.4, increases up 
to {q){z = 0) ~ 0.6, whereas d{q)/dr ~ 0. Dissipation 
enhances that process (Kazantzidis et al. 2004, Springel 
et al. 2004, Baffin 2005), specially in the innermost re¬ 
gion of the galaxy (r < O.Srrvir), although hydrodynami- 


Fig. 14.— Left column: Radius evolution. Top panel: Averaged 
Lz normalized to Lch obtained from the present distribution of 
debris as a function of stream age. Upper-middle panel: Averaged 
energy normalized to Ech as a function of stream age. In these 
panels open dots denote simulations in a halo with q = 0.8 = const., 
whereas filled dots denote debris in a halo with axis-ratio that 
increases linearly in time from q = 0.6 up to g = 0.8. In the 
bottom and middle-bottom panels we show the same quantities for 
satellites evolving in a halo with axis-ratio that decreases linearly 
in time from q = 1.0 down to g = 0.8 (filled dots) and a halo with 
g = 0.8 = const, (open dots). In all panels, full and dotted lines 
represent the progenitor’s quantities for a varying and a constant 
g, respectively. The initial and final mass of the satellite model are 
Ms{tQ) = IO^Mq and Ms{tf) = 5 x 10 ^Mq, respectively. Note 
that all simulations use an evolving halo. 


cal simulations appear to suffer from over-cooling, which 
over-estimates dissipation effects. 

In order to address whether the halo shape evolution 
can be uncovered from tidal streams, we have repeated 
our simulations in halos with a variable axis-ratio. For 
simplicity, we have assumed that q varies linearly with 
time from q{to) = 0.6 up to q{tf) = 0.8 (i.e., the halo 
becomes more spherical with time) and from q{to) = 1.0 
down to q{tf) = 0.8 (i.e, the halo becomes oblate with 
time). 

In Fig. in we show the radius evolution (left column) 
for both cases. Additionally, we plot the evolution of 
the model M1H2Q1 (i.e with q = 0.8 =const.) for com¬ 
parison. As expected from Fig. ^ increasing the value 
of q induces a reduction of the satellite’s peri and apo- 
galactica, and vice versa. 

In the right column we plot the present distributions of 
energy and the angular momentum of tidal debris as a 
function of stream age. Filled and open dots denote de¬ 
bris in a halo with varying and a constant axis-ratio, 
respectively, whereas the error bars represent the vari¬ 
ance around the mean value. Full and dotted lines show 
the look-back evolution of the progenitor’s energy and 
angular momentum. 

The results that we obtain are clearly equivalent to those 
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of H6.3I Although the orbit of the stream progenitor de¬ 
pends on the evolution of the halo potential, the present 
distribution of debris only reflects the present Galaxy 
potential . In particular, Fig. d shows the energy and 
angular momentum distributions are insensitive to the 
evolution of the halo’s shape if we force our models to 
end up with the same mass, position and velocity. 

8. EFFECTS OF DARK MATTER CLUMPS ON TIDAL 
STREAM PROPERTIES 

Having explored the effects of secular evolution of the 
smooth host potential, we now turn to the role of dark 
matter clumps in dispersing unbound (stream) particles. 
For simplicity, we assume that the progenitor’s orbit re¬ 
mains exclusively determined by the smooth component 
of the Galaxy potential, so that /sub lEo.lT^ only applies 
to unbound satellite particles. 

Note that in self-consistent simulations, the satellite’s 
orbit can also be altered through strong satellite-subhalo 
interactions, which would be easily detectable in the 
present stream’s energy-angular momentum distribution 
as discontinuities at different stream ages"^. Therefore, 
our approach minimizes the effect of dark matter clumps 
and only analyzes the heating induced by the halo’s 
lumpiness on tidal streams. 

We have carried out simulations with nine different sets 
of subhalos (see m drawing on the method in ^2.31 to 
construct merger trees and assuming an isotropic accre¬ 
tion of subhalos. Thus, for consistency, we consider only 
a spherical dark matter halo which evolves according to 
Eq. 0 and 0. 

For the parent satellites, we consider the orbits of 
the models M1H2Q2 and M2H2Q2. For each merger 
tree, we evolve these orbits and analyze the present- 
day E and distributions. Since the time required 
for the subhalo force calculation scales as where Ng 
is the total number of subhalos, and since in AGDM 
the number of subhalos in a given mass range goes as 
ns{M)dM ^ M'^dM, our study is limited to subha¬ 
los with masses above 1 O^M 0 (which implies a typical 
logio Ns i) 

By construction, the only effect of repeated encounters 
between tidal stream particles and subhalos is a progres¬ 
sive broadening of the stream distribution in the E, 
plane. In Fig. El we have plotted the ratio between the 
variance and the averaged value of the energy and the an¬ 
gular momentum as a function of the stream age for the 
models M1H2Q2 (left column) and M2H2Q2 (right col¬ 
umn) in the absence (dots) and presence (lines) of dark 
matter clumps. We find that: 

(i) As a result of encounters between subhalos and stream 
particles, the angular momentum dispersion shows a 
clear increase in the trailing tail (upper-left panel, dot¬ 
ted line) for age> 5 Gyr, whereas the angular momentum 
distribution of the leading tail (full line) has a dispersion 
that remains practically unaltered. That increase can be 
as large as a factor 1.75 for the model Ml, reducing to a 
factor 1.07 for the model M2. Remarkably, both models 
show that the stream heating can only be observed for 

^ We refer the reader to Pefiarrubia & Benson (2005) for a sta¬ 
tistical analysis of the effects and likelihood of clump—clump inter¬ 
actions in a Milky Way-like galaxy. These authors find that the 
present population of subhalos move on orbits that, statistically, 
have been barely altered by subhalo-subhalo interactions. 



Age (Gyr) Age (Gyr) 


Fig. 15. — Variance of the present energy (bottom panels) and 
angular momentum (upper panels) distributions as a function of 
stream age. Filled and open dots denote, respectively, the vari¬ 
ance of the leading and trailing tails in the absence of dark matter 
clumps. Solid and dotted lines show, respectively, the variance of 
the leading and trailing tails in galaxies with a subhalo population 
corresponding to that predicted by ACDM cosmology for a Milky 
Way-like galaxy. 

those parts older than 5 Gyr (i.e., t <9 Gyr). 

(ii) We do not observe any increase of as induced by 
subhalos. This results is in agreement with that found 
by Ibata et al. ( 2002 ), who showed that stream heating 
resulting from subhalo interactions is mostly reflected in 
the angular momentum distribution, whereas the energy 
distribution is barely altered. 

The angular momentum dispersion only shows a re¬ 
markable increase for age>5 Gyr. Following the results 
of Pefiarrubia & Benson (2005)- see also Zentner & Bul¬ 
lock (2003)- the accretion rate of substructures found a 
maximum at t = 4 Gyr, decreasing a factor 3 until t = 9 
Gyr. From t = 9 Gyr to the present, the accretion rate 
was approximately constant. 

This appears to indicate that the stream heating is cor¬ 
related with the accretion rate of subhalos: We observe 
heating only for age> 5 Gyr because the encounter rate 
between stream particles and dark matter clumps was ap¬ 
proximately a factor 3 larger for those particles stripped 
at t < 9 Gyr (age>5 Gyr) than for those that became 
unbound at t > 9 Gyr (age<5 Gyr). 

Finally, we must remark that the halo’s lumpiness will 
be more easily detected in the old (age> 5 Gyr) and cold 
(trailing tails of progenitors that were initially low mass) 
parts of tidal streams. 

9. DISCUSSION AND CONCLUSIONS 

In this paper we have analyzed what information can 
be extracted from stellar tidal streams on (i) the Milky 
Way’s halo shape, (ii) the halo’s secular evolution, (iii) 
the mass evolution of the stream’s progenitor and (iv) 
the presence of dark matter clumps in our Galaxy. We 
assumed as a boundary condition of this analysis that the 
present-day position, velocity and mass of the stream’s 
progenitor can be measured and are identical for all evo¬ 
lutionary scenarios. Under these conditions, we have ex- 
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plored whether the extended tidal debris reflect differ¬ 
ences arising from the items (i)-(iv) above. We have 
carried out this study for a Sagittarius-like dwarf galaxy, 
although our results are general and can be applied to 
other systems in the Milky Way. 

The main result is that tidal streams do not provide 
information on the adiabatic evolution of the Milky Way 
or, in other words, the properties of entire tidal streams 
only reflect the present Galaxy potential. Thus, ground- 
based observations already available for tidal streams 
(basically providing distances and radial velocities along 
the stream) and future satellite data covering the full 
phase-space (making possible studies in the E-L^ plane) 
can only constrain the present characteristics of the 
Milky Way potential. As a direct consequence. Galaxy 
evolution processes can be neglected when modeling tidal 
streams, which clarifies one of the main caveats in cur¬ 
rent N-body simulations and confirms that tidal stream 
models computed under that hypothesis are appropriate 
for tracing the distribution of dark matter around our 
Galaxy. 

In contrast, we And that tidal streams are indeed sen¬ 
sitive to the present properties of halo shape. In par¬ 
ticular, we confirm that measuring the precession rate is 
a fairly powerful method to constrain the halo flatten¬ 
ing of the gravitational potential. This may not require 
measurements of proper motions. 

We have shown that the study of tidal streams in 
the E-Lz plane provides information on the progenitor’s 
mass loss since the time of accretion. The energy-angular 
momentum distribution of stream particles has an av¬ 
erage value that only reflects the present position and 
velocity of the progenitor system. In contrast, the E- 
Lz variance about that mean increases with the initial 
satellite mass, thus, making it possible to determine the 
mass loss fraction directly from the present E-Lz dis¬ 
tribution. Furthermore, since a secular drift in energy 
and angular momentum is induced by dynamical friction 
(a drag acceleration that scales in proportion to Mg) one 
could, in principle, reconstruct the mass loss curve Ms{t) 
from the age of different stream pieces (if the age can be 
labeled e.g. by metalicity or by theoretical modeling), 
which ultimately depends on the initial mass profile of 
the satellite galaxy (see Zhao 2004). 

We have analyzed the effects of dark matter clumps 
on tidal streams in Section We simplify the problem 
by assuming that dark matter clumps do not alter the 
progenitor’s orbit, but only the orbits of stream parti¬ 
cles. That approach establishes, therefore, a minimum 
impact of subhalos on tidal stream properties (we note 
that a sharp change of the progenitor’s Ep^ Lz,p induced 
by a collision with a subhalo at t = tc would be re¬ 
flected as a discontinuity in the averaged stream’s E, Lz 
at age=tf—tc). We have confirmed that dark matter sub¬ 
halos induce only very modest stream “heating” by in¬ 
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creasing of the angular momentum dispersion in the old¬ 
est (age>5 Gyr) and coldest (trailing tail) stream parts. 
This raises the question of whether one can constrain 
the halo lumpiness either with current ground-based 
techniques or with future astrometric satellite missions 
(GAIA, SIM). 

At present, the detection of the oldest parts (i.e those 
stars that became unbound first) of a tidal stream with 
state-of-the-art ground-based surveys is challenging. Al¬ 
though the trailing tail maintains a coherent structure 
and should be easier to detect as spatial over-densities 
(see Sec. Iti.4ll . its surface brightness decreases consider¬ 
ably with time, which reduces the possibility of detec¬ 
tion above the Galactic field contamination in large field- 
of-view color-magnitude diagrams like those provided 
by SDSS. Furthermore, it is also expected that tidal 
streams are composed by old, metal-poor stellar pop¬ 
ulation. Therefore, some valuable techniques for tracing 
tidal streams with all-sky surveys (like 2MASS) cannot 
be applied to detect the oldest stream pieces since these 
are barely sensitive to metal-poor stars expected in old 
parts of tidal streams. Also surveys of tidal streams using 
M-giant stars (Majewski et al. 2004) are limited to the 
youngest stream pieces (material unbound only 1-2 Gyr 
ago; LJM), while the oldest wraps still remain hidden in 
the Galactic halo. Current theoretical models of the two 
largest, brightest streams in the Milky Way (Sgr: LJM; 
Monoceros: Peharrubia et al. 2005) indicate that there 
is no detection of tidal debris that became unbound more 
than 2-3 Gyr ago. Finally, stream tails are not located 
on the progenitor’s orbital plane (see Sec. lOl . which 
further complicates tagging of debris as part of a known 
stellar system and needs of accurate stream models in 
order to determine a possible common origin and to es¬ 
timate the stream’s age. 

In the future, the most powerful method to search for an¬ 
cient debris in the halo should come with the next gener¬ 
ation of astrometric satellites, which will permit analysis 
of tidal streams in the constant-of-motion space, thus, 
providing the most straightforward way to identify tidal 
debris independently of the stream age. However, it is 
important to remark that observational errors may in¬ 
troduce strong limitations, as indicated by Brown et al. 
(2005). Possibilities of identifying satellite remnants in 
the Milky Way halo are considerably reduced after tak¬ 
ing into account the observational errors expected for 
the GAIA catalogue and the large number of background 
stars. 
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